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In an ideal hydrodynamic model, with an equation of state where the confinement-deconfinement 
transition is a cross-over at T^o = IQ&MeV , we have simulated v^=200 GeV Au+Au collisions. 
Simultaneous description of the experimental charged particle's pr spectra and elliptic flow require 
that in central (0-10%) Au-f Au collisions, initial energy density scales with the binary collision 
number density. In less central collisions, experimental data demand scaling with the participant 
density. Simulation studies also indicate that in central collisions viscous effects are minimal. 
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One of the important findings in Au-f Au collisions at 
RHIC is that the centrality dependence of particle mul- 
tiphcity can be understood in a simple geometric model. 
For example in it was assumed that a fixed fraction 
X of multiplicity per unit (pseudo) rapidity Upp 



drj 



in pp collisions is due to hard 'processes' and the rest 
(1 — x) is due to soft processes. Assuming that hard pro- 
cesses scales with the binary collision numbers [Ncoii) 
and the soft processes scales with the participant num- 
bers (Npart), the pseudo-rapidity density in nucleus- 
nucleus collisions is then parameterised as, 



dN 
drj 



part 



(1) 



Using Glauber model to compute Npart and Ncoii, 
PHOBOS data 3 O'^ charged particle multiplicity, in 
Au-f Au collisions were fitted with hard scattering frac- 
tion X = 0.05 ± 0.03 at yi=56 GeV and x = 0.09 ± 0.03 
at Vs=130 GeV Q. More recently, PHOBOS collabo- 
ration studied the geometric scaling of (pseudo) rapidity 
density in ^5=19. 6 and 200 GeV Au-f Au cohisions 
For both the collision energies, hard scattering fraction is 
approximately constant, x = 0.13±0.01(stai)±0.05(sys). 

Geometric scaling of Au-f Au collisions as in EqlU has 
been widely used in hydrodynamic model calculations [3- 
[6|. Hydrodynamic models require initial energy density 
configuration. Following Eq[Tl in an impact parameter b 
Au-f Au collision, initial energy density in the transverse 
plane can be parameterised as. 



e(b, x, y) = eo[(l-x)Npart{^,x,y)+xNcou{h,x,y)], (2) 

where NpartO^TX,y) and Ncou{h,x,y) are the transverse 
density distribution for the participant pairs and the col- 
lision number respectively, x is the fraction of hard scat- 
tering. Parameterisation Eql5]is generally called Glauber 
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model initialisation. Eq in Eql2]is the central energy den- 
sity in b = collision. It is generally fitted to experimen- 
tal data, e.g. multiplicity, pr spectra etc. Ideal hydro- 
dynamics, with Glauber model initialisation, with hard 
scattering fraction x — 0.25 explains a variety of experi- 
mental data, e.g. identified particle's multiplicity, mean 
momentum, pT spectra, elliptic flow etc |5|. Ideal hy- 
drodynamics description to the experimental data how- 
ever deteriorates beyond pr «1.5 GeV, presumably due 
to increasing role of dissipative effects at large pt- De- 
scription also deteriorates in very peripheral collisions. 
Glauber model initialisation of the energy density, with 
X = 0.13, also give reasonable description to the exper- 
imental data [6|. One inconsistency however remained. 
Glauber model initialisation with hard scattering frac- 
tion X = 0.25 or a; = 0.13, under predict experimental 
elliptic flow in very central, e.g. 0-10% collisions. 

Elliptic flow is a key observable in establishing that 
in ^=200 GeV Au-fAu collisions at RHIC, lattice 
QCD predicted Quark-Gluon-Plasma (QGP) is pro- 
duced. QGP is a thermal system. Finite elliptic flow in 
Au-fAu collisions and the fact that hydrodynamic mod- 
els do explain the flow are generally cited as proof of QGP 
production in Au+Au collisions. It is then important to 
understand why Glauber model of initial condition un- 
derestimate elliptic flow in very central collisions. Neglect 
of dissipative effects can not be the reason. Inclusion of 
dissipative effects only reduces elliptic flow. 

To understand the relation between elliptic flow in cen- 
tral Au-fAu collisions and the geometric scaling of initial 
energy density as in EqUJ we have simulated -^5=200 
GeV Au-fAu collisions with Glauber model initial condi- 
tion at two extreme limits of the hard scattering fraction, 
X = and a; = 1 . Two limits corresponds to very different 
collision dynamics, for x = 0, initial energy density scales 
with participant density implying that Au-fAu collisions 
are completely dominated by the soft processes, while 
for X = 1, Au-fAu collisions are completely dominated 
by the hard processes and the initial energy density scales 
with the density of binary collision numbers. Actual sce- 
nario may be in between the two extreme limits. 

Details of a hydrodynamic model can be found in Q . 
We assume that in Au-fAu collisions, a 'baryonless', 
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'ideal' QGP fluid is produced. Space-time evolution of 
the fluid is obtained by solving the energy-momentum 
conservation equation, 

with the code 'AZHYDRO-KOLKATA' in (r = 

y/W^^z^, x,y,7y = ^Inf^f) co-ordinate system. Lon- 
gitudinal boost-invariance is assumed. Eq|3] is closed 
with an equation of state (EOS), p = p(e). We use the 
recently constructed lattice-|-HRG EOS where the 
confinement-deconfinement transition is a cross-over at 
Tco^ige MeV. The deconfined part of the EOS corre- 
sponds to recent lattice simulations while the con- 
fined part corresponds to that of a hadronic resonance 
gas with all the resonances with mass irises <2.5 GeV. 
In 11] it was shown that the lattice based EOS with 
confinement-deconfinement cross-over transition, reason- 
ably well explain the centrality dependence of </> mesons 
multiplicity {dN'^ / dy) , mean pr {{Px)) ^-nd integrated 
elliptic flow (^2 )• Recently, it was shown that hydrody- 
namical evolution of QGP fluid, with the lattice-|-HRG 
EOS reasonably well explains the transverse momentum 
spectra of tt, K and </> mesons j^lSj]. Proton spectra how- 
ever are under predicted. 

Solution of Eq|3] require initial conditions, e.g. the ini- 
tial or the thermalisation time beyond which hydrody- 
namics is applicable, initial energy density, fluid velocity 
etc. We assume that the fluid is thermalised in the time 
scale Ti=0.6 fm 0. At T,=0.6 fm, energy density in the 
transverse plane is distributed as in Eq|31 with hard scat- 
tering fraction (i) cc = or (ii) x = 1. Irrespective of the 
hard scattering fraction, initial fluid velocity is assumed 
to be zero, Vx{x,y) = Vy{x,y) = 0. Hydrodynamic mod- 
els also require a freeze-out prescription to convert the 
information about the fluid energy density and velocity 
to particle spectra. We assume that the fluid undergoes 
kinetic freeze-out at temperature rF=150 MeV (for bary- 
onless fluid chemical potential is zero throughout the evo- 
lution). 

With hard scattering fraction fixed, either x — or 
X = I, central energy density Eq is the only parame- 
ter left in the model. We fix eo by fitting the PHENIX 
data on charged particles pr spectra in 0-10% Au-f Au 
collisions [l3]. For x = 0, best fit to the 0-10% data is 
obtained with eo=36.1 GeV/fw?. PHENIX data require 
~30% higher central energy density, eo=48 GeV/ fm?, if 
Au-|-Au collision is completely dominated by the hard 
processes {x = 1). The solid and the dashed lines in 
FigH^ are fit to the data with hard scattering fraction 
X = 1 and X = Q respectively. Data are well fitted. 
Charged particles pr spectra in central Au-|-Au colli- 
sion is insensitive to the hard scattering fraction x in 
the Glauber model of initial condition. 

With initial energy density fixed, we can predict for 
the pt spectra in all the other collision centralities. In 
Figlll in panels (b-f), model predictions for the charged 
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FIG. 1: In panels (a-f) PHENIX data 14] on the charged 
particle transverse momentum spectra in 0-60% Au+Au col- 
lisions are compared with hydrodynamic model predictions. 
The dashed lines are the predictions when initial energy den- 
sity scales with participant density (a; = 0). The solid line are 
predictions when energy density scales with density distribu- 
tion of binary collision numbers (x = 1). 



particles pr spectra in 10-20%, 20-30%, 30-40%, 40-50% 
and 50-60% Au-|-Au collisions are compared with the 
PHENIX data [13] • When central energy density is fixed 
to reproduce charged particles pt spectra in 0-10% colli- 
sions, Glauber model initialisation with x = 1, also give 
reasonable description to the spectra in 10-20% and 20- 
30% centrality collisions. But in more peripheral colli- 
sions, PHENIX data are under predicted. Initialisation 
with hard scattering fraction x = however continue 
to explain data till very peripheral collisions, though 
description to the data deteriorates at larger pT or in 
more peripheral collisions. It appear that if in Au-|~Au 
collisions, energy density scales with participant density 
(x = 0), charged particles px spectra in 0-60% Au-|-Au 
collisions are reasonably well explained. On the contrary, 
if energy density scales with binary collision number den- 
sity (x = 1), charged particles p^ spectra, only in central 
(0-30%) collisions is explained. 

Let us now study centrality dependence of simulated 
elliptic fiow. In a hydrodynamic model, elliptic fiow de- 
pends on the initial spatial eccentricity. Ex = Ifnp^, (...) 
denotes energy density weighted averaging. In Fig|51 cen- 
trality dependence of Ex , in the two extreme limits x = 
(filled circles) and x = 1 (filled squares) are shown. Ex 
is more if Au-|-Au collisions are dominated by the hard 
processes (x = 1) rather than the soft processes (x = 0). 
Glauber model initialisation of energy density with hard 
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FIG. 2: Initial eccentricity of the collision zone as a function 
of participant numbers. The filled circles and squares are 
eccentricity with hard scattering fraction and 1 respectively. 



scattering fraction x = 1 will generate more elliptic flow 
than the initialisation with x — 0. 

In FiglSl we have compared the simulated V2 with 
PHENIX measurements [TB]. PHENIX collaboration ob- 
tained V2 from two independent analysis, (i) event plane 
method from two independent subdetectors, V2{BBC} 
and V2{ZDC — SMD} and (ii) two particle cumulant 
W2{2}. i'2{2} from two particle cumulant and V2{BBC} 
or V2{ZDC — BBC} from event plane methods agree 
within the systematic error. All the measurements upto 
Pt = 3 GeV are shown in FiglS] The solid and dashed 
lines in FigE] corresponds to Glauber model initial con- 
ditions with X — 1 and x = respectively. As expected 
from the eccentricity study (Figl2]), in all the collision 
centralities, Glauber model initialisation with x = 1 gen- 
erate more flow that the initialisation with a; = 0. Unlike 
the pt spectra in central Au-f-Au collisions, which do not 
distinguish between the initial conditions with x = 1 and 
a; = 0, elliptic flow, being a more sensitive observable, 
can distinguishes between them. It is very interesting to 
note that Glauber model initialisation with hard scat- 
tering fraction a; = I well explain the PHENIX data on 
elliptic flow in 0-10% Au-f Au collisions. However, in all 
the other collision centralities elliptic flow is over pre- 
dicted. For example, at pr «1.5 GeV, simulated flow 
with X = 1 over predict experiments by ^20%, 25%, 
35%, 45% and 60% in 10-20%, 20-30%, 30-40%, 40-50% 
and 50-60% Au-|-Au collisions. At larger pT flow is even 
more over predicted. Glauber model initialisation with 
hard scattering fraction x = predict less flow and el- 
liptic flow in 0-10% Au-f Au collisions is underpredicted 
e.g., at Pt ~1.5 GeV, it under predict the experiment 
by ^35%. In all the other collision centrality agreement 
with data is comparatively better. In 10-20% and 20-30% 
Au-|-Au collisions, Glauber initialisation with a; = give 
very good description of the data up to pr ~1.5 GeV. In 
more peripheral collisions, flow is over predicted. Even 
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FIG. 3: The filled circles, squares and triangles are PHENIX 
measurements 15] for elliptic flow in 0-60% Au-|-Au collisions. 
The solid and dashed lines are elliptic flow in hydrodynamic 
simulations with hard scattering fraction x=l and x=0 re- 
spectively. 



then agreement with data is better than that obtained 
with hard scattering fraction a: = 1. For example at 
Pt wl.5 GeV, simulated flow over predict PHENIX data 
by - 15%, 25% and 50% in 30-40%, 40-50% and 50-60% 
collisions. 

Considering that flnite elliptic flow is an essential ob- 
servable to establish that QGP is produced in -^5=200 
GeV Au-j-Au collisions, it is imperative that hydrody- 
namical models, with QGP as initial fluid, explains the 
flow in 0-10% Au-|-Au collisions in addition to the pT 
spectra. Present analysis indicate that in 0-10% Au-l-Au 
collisions, simultaneous description of pt spectra and el- 
liptic flow require hard scattering fraction a; = 1 in the 
Glauber model of initial condition. In less central col- 
lisions however, simultaneous description of pt spectra 
and elliptic flow are best obtained with hard scattering 
fraction a; = 0. The result implies that geometric scal- 
ing of Au-|-Au collisions changes with collision centrality. 
In a central collision, energy density scales with binary 
collision number density while in a less central collision, 
energy density scales with participant density. Arguably, 
transition from binary collision number scaling to par- 
ticipant scaling can not be as sharp as conjectured here. 
More detailed analysis is required to flnd the width and 
exact location of the transition. 

Such a transition may have implications for the hydro- 
dynamical analysis also. Note that for binary collision 
number scaling, fluid has to be initialised at higher en- 
ergy density (£o=48 GeV/ fm? for a; = 1 and eo=36.1 
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GeV/ fm^ for x = 0). We have assumed similar thermal- 
isation time ri=0.6 fm for both the scahng conditions. 
Since thermalsation time scale is expected to be inversely 
proportional to the density, it is likely that the fluid in 
central Au+Au collisions will thcrmalise in a lesser time 
scale that the fluid in less central collisions. The re- 
sult may also have implication on the dynamics of the 
pre- equilibrium stage. Hydrodynamic models assume lo- 
cal thermalisation. The fluid produced in Au-I-Au col- 
lisions evolves through a pre-equilibrium stage to equi- 
libration. At present, we have limited knowledge about 
the pre-equilibrium stage. Present results suggests that 
in 0-10% Au-|-Au collisions, pre-equilibrium stage is dom- 
inated by binary collisions, but in a less central collision, 
pre-equilibrium stage is dominated by the 'wounded' nu- 
cleons. 

In the present analysis, the effect of viscosity is ne- 
glected. QGP viscosity is a contentious issue. String 
theory based models (ADS/CFT), which are unrelated 
to QCD, give a lower bound on viscosity of any mat- 
ter, viscosity to entropy ratio, f]/s > l/in [l6l|. Per- 
turbative estimate of the viscosity is also uncertain to a 
great extent, rj/s 0.0-1.0 |17H19| . Several authors have 
extracted viscosity directly from the experimental data 
[ni [13413 • Depending on the model, data analysed etc. 
the extracted viscosity vary over a large range and one 
can possibly give an upper bound, rj/s < 0.5 i27||. Ef- 
fect of viscosity is to enhance particle production mainly 
at large px, and also to reduce elliptic flow. One observes 
from FiglHandO that in 0-10% Au-|-Au collisions, if ini- 
tial energy density scales with collision number density, 
ideal fluid dynamics hardly leave any scope for viscous 
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